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Abstract 



Vacuum polarization of massive scalar fields in a thermal state at arbi- 
trary temperature is studied near the horizon of a Reissner-Nordstrom black 
hole. We derived an analytic form of (eft 2 ) approximately in the large mass 
limit near the black hole horizon. We uses the zeroth order WKB approxima- 
tion and power series expansion near the horizon for the Euclideanized mode 
function. Our formula for the vacuum polarization shows regular behavior 
on the horizon if the temperature of the scalar field is equal to the Hawking 
temperature of the black hole. The finite part of the vacuum polarization 
agrees with the result of the DeWitt-Schwinger approximation up to 0(m~ 4 ) 
which is the next leading order of the expansion. 
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I. INTRODUCTION 



Since Hawking's discovery of black holes radiation M, it has been an important subject 
to investigate quantum effects in black hole spacetimes and a great deal of effort has been 
done. The expectation values of the regularized vacuum polarization ((f) 2 ) and the stress 
energy tensor (T^) are important quantities to understand quantum effects in black hole 
spacetimes. They represent the effect of quantum fluctuations and back reaction on black 
hole spacetimes and play an important role in the context of semi-classical theory of quantum 
gravity. 

For a scalar field in the large mass limit, of which Compton length is much smaller 
than the characteristic curvature scale of background geometry, we can evaluate regularized 
((f) 2 ) and (Tfj_ u ) by using the DeWitt-Schwinger approximation B-0]. This is an adiabatic 
expansion in terms of inverse power series of the mass of scalar field. One can evaluate 
divergent and finite parts of ((f> 2 ) and (T^ u ) up to arbitrary order. Frolov and Zelnikov |||| 
suggest that the finite parts of ((f> 2 ) and (T^ v ) calculated by the DeWitt-Schwinger expansion 
coincides with the value calculated by assuming the Hartle-Hawking state for the scalar 
field, although the DeWitt-Schwinger expansion does not impose any condition on the state 
of the scalar field. The justification is confirmed numerically near the horizon |§|§. 

The Hartle-Hawking state is defined by requiring that the Euclidean Green's function is 
regular on the horizon 0. Its physical meaning is that a black hole is placed in a cavity 
and is equilibrium with thermal black-body radiation. The thermal state is defined by the 
Euclidean Green's function by imposing periodicity in the imaginary time coordinate with 
a period of the inverse of the temperature. On the order hand, in the DeWitt-Schwinger 
expansion, the condition which characterizes the thermal state is not imposed. So it is 
not obvious whether the DeWitt-Schwinger expansion yields the same finite part of the 
Hartle-Hawking state near the horizon. 

In this paper, we consider scalar fields in static spherically symmetric background space- 
times. By assuming that the mass of scalar fields is large, we calculate ((f> 2 ) in the thermal 
state with arbitrary temperature using point splitting method near the horizon. A similar 
calculation is already done by Anderson et al. but their expression of the vacuum po- 
larization diverges on the horizon which contradicts with their numerical calculation. Our 
purpose is to derive the vacuum polarization which is regular on the horizon, and resolve 
why the finite part of ((f) 2 ) in the Hartle-Hawking state near the horizon can be approximated 
by the DeWitt-Schwinger expansion which does not assume any condition on the thermal 
state. 

The plan of this paper is as follows. In Sec. II, we review the method to calculate ((f) 2 ) 
in a thermal state using point splitting regularization. In Sec. Ill, we present our method 
to calculate ((f) 2 ) in the large mass limit near the black hole horizon. Sec. IV is devoted to 
conclusion and discussion. Throughout the paper, we use units such that h = c = G = ks = 



1. Our sign conventions are those of Misner, Thorne and Wheeler [10 



II. DERIVATION OF VACUUM POLARIZATION (<p 2 ) IN A THERMAL STATE 

We consider a minimally coupled scalar field with mass m in a thermal state at arbitrary 
temperature T in a Reissner-Nordstrom spacetime. The metric in Euclidean section r = —it 



2 



is given by 



ds 2 = fdr 2 + jdr 2 + r 2 dn 2 , 



(1) 



where / = (r — r + )(r — r_)/r 2 , r + and r_ are the location of inner and outer horizon, 
respectively. The surface gravity of a black hole is given by k = (r + — r_)/(2rl). The 
vacuum polarization (0 2 ) can be computed from the Euclidean Green's function by noting 
the identity 



2 )unreg = T:G (1 \x,x) = lG F (x,x) = G E (x,x), 



(2) 



where G^ 1 ' is the Hadamard Green's function, Gp is the Feynman Green's function, and 
Ge is the Euclidean Green's function, respectively. Ge is divergent at coincident limit and 
must be regularized. In this paper, the covariant point-splitting method is employed for the 
regularization of ultraviolet divergences. We start from an expression for G E (x,x') when 
the points x and x' are separated. Next we prepare an appropriate counterterm to subtract 
the divergence, and then take coincident limit. This procedure is shown schematically as 



2 ) reg = Jim [G E (x, x') - G DS (x, x')} 



(3) 



where Gds( 

point splitting counterterm needed to regularize G E (x,x'). The Eu- 
clidean Green's function obeys the equation 

S 4 (x, x') 



[D-m 2 ] G E (x,x') 
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where m is the mass of the scalar field and d'Alembertian □ evaluated using the Euclidean 
metric ([I]). The point splitting counter term needed to renormalize Ge{x,x') in arbitrary 
spacetime is 



G DS {x,x') 



1 , m 2 + (£-±)i2 
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where ^ is coupling constant to the scalar curvature R, a is the one half the square of 
geodesic distance between the points x and x' and 7 is Euler's constant. 

For simplicity of calculation, the quantity ((f) 2 ) is evaluated by timelike point separation 
t' 7^ t. By assuming that the Green's function is periodic in e = r — r' with period 1/T, the 
expression for the vacuum polarization in the thermal state at temperature T is given by 

(cj) 2 ) T = lim [G E (x, r; x, r') - G DS (x, r; x, r )] 
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where 
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The term l/(r\Z/) is subtracted to cancel the superficial divergence which comes from the 
choice of timelike point-splitting. p n i and q n i are independent solutions of the mode equation 
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and satisfy the Wronskian condition 

dq n i 
dr 
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Anderson et al. [|8||| calculated the second and the third term of 
the second order WKB approximation . Their result is 

m 2 f 
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]) analytically using 
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The expression (|T0|) was obtained by taking into account only the contribution from n > 1 
mode. In the case of massless scalar fields, excluding the lowest frequency mode n = 
corresponds to imposing an infrared cutoff. The expression ( |T0"D gives good approximation for 
massless fields both near and far from the horizon in Schwarzschild and Reissner-Nordstrom 
spacetimes. This was confirmed by their numerical calculation pi. 

When the result (|T0|) is applied to the massive scalar field, however, two unpleasant 
features appear. Firstly, the regularized expression always has logarithmic divergence on the 
horizon / = even if the scalar field is in the Hartle-Hawking state. It is not consistent with 
their numerical work which shows regular behavior of {4> 2 )t on the horizon. Secondly, the 
finite terms of 0(m 2 ) and 0(m°) are left in the regularized expression. This disagrees with 
the result of the DeWitt-Schwinger approximation by Frolov [|] which does not contain such 
finite terms. Therefore the expression fllCf ) cannot yield good approximation for the massive 
scalar field. In the next section, we present the method to improve their approximations 
and resolve these problems. 



III. VACUUM POLARIZATION NEAR THE HORIZON IN THE LARGE MASS 

LIMIT 

In this section, we present the method to resolve the unwanted behavior arose in the 
expression of (4> 2 )t near the horizon / ~ 0. We assume that the mass m of the scalar field 
is large: 

m>— . (11) 
We evaluate n > 1 and n = contribution to (<P 2 )t separately. 
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A. n > 1 contribution 

We can use the WKB approximation for n > 1 mode. The zeroth order WKB solution 
is given by 
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The WKB approximation is correct near the horizon. Using this solution, we have 

rp oo 

{4> 2 )n>\ = — ^cos(nfte)/ n , 



n=l 



(13) 



where 




(14) 
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J n can be evaluated using the Plana sum formula |8||| under the assumption / ~ 0, a r 
Knjf: 
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Therefore, for n > 1, 



1 2 / 2 r 2 n 2 n 2 
~+r fm +— — 



1/2 



+ 



/ 



rftra 



ftn m 



/ 2ftn 



rn 3 r 2 nn 



(16) 



5 



By taking n sum, 0(f /(rn 3 ), \fj /(r 2 K,n)) terms give finite values which becomes zero on 
the horizon and we can omit these terms. Using the formula 



cos(nKe) 1 2 2 2 , 
^-^ = --ln( K 2 e 2 ) + 0(e 2 ) 

n=l - 

1 K 2 



(17) 



k Tin cos(nKe) = — \- 0(e 



n=l 



we have 
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Using the DeWitt-Schwinger counter term ([/]), the regularized expression of the expectation 
value becomes 



m 



( m 2 f \ m 2 T 2 f' 2 f" f 



(0 )n>l - (0 }ds - In j + ^ + — - - - (19) 

This gives the same expression of Anderson et al. ||, who used the second order WKB 
solution. Near the black hole horizon, it is sufficient to use the zeroth order WKB solution 
to reproduce the result of the second order WKB approximation. On the horizon / = 0, 
this expression diverge and we must examine the contribution from n = mode. 



B. n = contribution 

We must recall that the expression (p~Q|) does not include n = mode but there is no 
reason to exclude this contribution in the case of the massive scalar field. So we must 
investigate the contribution of n = mode. Since the WKB method breaks down for n = 
mode near the horizon, we cannot apply the approximation used by Anderson et al. || . We 
solve n = mode function by the following method. 

As we assume that the mass of the scalar field is large, a dimensionless constant 

r+ ~ r - (20) 



4r + [m 2 r 2 + + 1(1 + 1)] 
becomes smaller than unity for all value of I: 

e < 1 (for all I) (21) 
We rescale the radial coordinate r as follows: 



x = J r —^. (22) 
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Assuming that x is 0(1) quantity, this rescaling of the radial coordinate means we are 
concentrate on the region (r — r + ) ~ 0(e), which is the region near the black hole horizon. 
Using x and e, the radial equation ([8]) for n = mode can be written as 



d 2 Xoi 
dx 2 



+ 



- + F(x,e) 
x 



0/ 



dx 



[l + G(x,e)} X oi = 0, 



(23) 



where 



F(x,e) 



G{x,e) 



2ex 



2nr + + ex 2 

16e 
2kt + + ex 2 



[nr + + 2m 2 r 2 + (2e 2 x 2 + 3e 



(24) 
(25) 



We can express the approximate solution of the mode equation ( |23| ) by the power series 
expansion with respect to a small parameter e: 



Xoi (x,v) = Xm {x) + e Xm ' (x) + e 2 Xoi' (z) + 



.2^(2), 



(26) 
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(27) 



(for n > 1) 



(28) 



We can obtain solutions of these equations by using modified Bessel functions. For example, 
the lowest order solution is given by 



(0) 
PlO 



Io(x) 



(0) 

9io 



-Kn(x). 



(29) 



Substituting the expansion (|26|) to 
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E 

1=0 



(21 + i)pi q m - 



(30) 



and taking summation about I using the Plana sum formula, (0 2 ) n=o is written as the power 
series expansion with respect to (mr + ) _1 . We calculated (0 2 ) n=o up to O ((mr + )" 4 ), which 
is accomplished by the e expansion of the mode function up to 0(e 3 ). By taking into the 
account of the contribution of n = mode, we finally get the form of the vacuum polarization 
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near the horizon as follows: 
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The divergences in (0 2 )t appears only in the terms of 0(m 2 ) and 0(m°). Taking / — > 
limit, 0(m 2 ) term diverge as ln(r — r+) and 0(m°) as (r — r+)~ 1 . These divergence disappear 
completely if the temperature of the scalar field equals T = k/2tt = Th, where T# is the 
Hawking temperature. Specifying T = X# is equivalent to require the regurality of (</> 2 ) on 
the horizon. In this case, the black hole is equilibrium with the thermal scalar field with 
temperature T = Tjj. This is nothing but a Hartle-Hawking state. 

We can evaluate a finite part of (0 2 )t i n the Hartle-Hawking state on the horizon. The 
leading order of {4> 2 )t at T = T# is 0{m~ 2 ) because not only divergent but also finite terms 
contained in 0(m 2 ) and 0(m°) terms vanish. The result is given by 



T=T„ 



+ 



7207r 2 m 2 r\ 
1 1 
25207r 2 m 4 r! 



6r_ 



+ 



4r z 



•10 



r_ 
39— 



2 3 

52^ + 23^ 



O (m" 6 ) . 



(32) 



This expression agrees with the result of the DeWitt-Schwinger expansion 



IV. CONCLUSION AND DISCUSSION 



We calculated (0 2 )t of a massive scalar field near the horizon of a Reissner- Nordstrom 
black hole using point splitting method. Our results are as follows: (i)all divergent and finite 
terms of 0(m 2 ) and O(m ) in (0 2 )t are canceled to be zero if the temperature of the scalar 
field equals to the Hawking temperature Th of a black hole, (ii)at this temperature, the 
leading order of the regularized (4> 2 )t becomes 0(m~ 2 ) and its value on the horizon agrees 
with that of the DeWitt-Schwinger expansion up to 0(m -4 ). (i) and (ii) indicate that regu- 
larized (4> 2 )t in the Hartle-Hawking state near the black hole horizon is well approximated 
by the DeWitt-Schwinger expansion. This confirms the results of [H>§|- 

It is not trivial that the finite part of (4> 2 )t i n the thermal state at T = Th is reproduced 
by the DeWitt-Schwinger expansion which does not specify any thermal state of the scalar 
field. It was shown that the leading order of the DeWitt-Schwinger expansion reproduces the 
value (0 2 )t=o of scalar field with large mass in general static spherical symmetric spacetime 
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TT| . This result is natural because the scalar field with T = is in a vacuum state and the 
DeWitt-Schwinger expansion does not incorporate information on a thermal state. 

Discrepancy between (0 2 )t in the Hartle-Hawking state and that in he DeWitt-Schwinger 
expansion will appear far from the horizon, because the thermal effect dominates the gravi- 
tational effect in the asymptotic flat region. If the thermal effect is dominant or comparable 
to the contribution from curvature effect near the horizon, there would also be discrepancy 
between (0 2 )t in the Hartle-Hawking state and that of the DeWitt-Schwinger expansion. 
However this is not seen from our results. This implies that the finite part of (4> 2 )t in 
the Hartle-Hawking state near the black hole horizon is dominated by the contribution from 
curvature effect. Why the contribution from thermal effects to finite parts can be negligible? 
The characteristic scale for thermal state is the inverse of the temperature 1/Th ~ M which 
is a scale length of thermal fluctuation. On the other hand, the Compton length of the 
scalar field is 1/m. The contribution of the Hawking radiation to (0 2 ) is oc Tfj/m? and that 
of thermal excitation is oc exp(— m/Tjj). For m 3> the effect of thermal excitation is 

exponentially small. Furthermore, near the black hole horizon, almost all part of thermally 
excited particle is absorbed by the black hole. This is the reason of the DeWitt-Schwinger 
approximation gives the same result of Hartle-Hawking state on the horizon. 
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